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RANDOM RESPONSE OF RECTANGULAR PANELS TO 


THE PRESSURE FIELD BENEATH A TURBULENT 
BOUNDARY LAYER IN SUBSONIC FLOWS 

Wei J. Chyu and M. K. Au-Yang 
Ames Research Center 


SUMMARY 


The response of a rectangular panel under the excitation of a turbulent boundary layer 
with a zero longitudinal mean pressure gradient in a subsonic flow was studied in detail. The 
method of normal mode was used together with the technique of spectral analysis. Simply 
supported and clamped edge boundary conditions of a panel were assumed. The response of 
the panel has been expressed in terms of displacement power spectral density. The mode number 
and frequency in the response computation were extended to 7 and 3000 Hz, respectively. 

The results of response computation assuming a clamped edge condition compare favorably 
with existing experimental data with a clamped edge condition. The simply supported edge 
condition, however, overestimates the response. The effect of flow velocity on the response is 
to increase the higher frequency components of the power spectral density per unit excitation. 

Charts of structural acceptances, which provide a framework for estimating the response of 
other rectangular panels, are presented as functions of dimensionless frequency and boundary- 
layer displacement thickness. Physical significances of the structural acceptances are also discussed, 
as regards the coincidence of the pressure and flexural waves and the probability of mode transition. 


INTRODUCTION 


The pressure fluctuations in regions of turbulent attached and separated boundary layers 
and shock waves adjacent to the surface of aerospace vehicles cause structural vibrations throughout 
atmospheric flight. The study of these vibrations is important in determining stress, fatigue, life 
of structures, and noise transmission into the interior of the vehicle. Unfortunately, the analysis 
of this type of vibration is complicated by the inherent random characteristics of the excitation 
pressure fields, and the difficulty of analytically describing the vibration of a realistic structure. 
For these reasons, early investigators of this problem considered only a hypothetical flow field 
and made the simplifying assumption that the structure, almost invariably either a beam or a 
rectangular panel, was infinitely large (refs. 1, 2, and 3). This assumption leads to a solution in 
terms of the mean square displacement of the panel as a whole, but not as a function of location 
on the panel. 



In recent investigations the response of finite-size rectangular panels to the excitation of an 
attached turbulent boundary layer has been considered, but the panels were assumed to be simply 
supported (refs. 4 to 6). This assumption simplifies the algebra tremendously, permitting solutions 
to be expressed in closed forms. Although these theoretical results agree better with experiment 
than those obtained with the infinite panel assumption, they tend to overestimate the response 
of a realistic panel. Previous analysis also fails to predict the wave matching accurately between 
the flexural wave of a panel and the pressure wave. The estimation of this wave matching is 
important, particularly when the matching occurs at one of the resonant frequencies of the 
structure, and causes a large structural response. With the advent of modern high speed computers, 
a theoretical analysis of a finite rectangular panel with clamped edges under the excitation of a 
turbulent boundary layer is now feasible. Analytical integration in closed forms is not necessary, 
since numerical integration can be carried out with no algebraic simplification of the integrand. 
This digital computer oriented approach has the following additional advantages over an analytical 
approach: (1) The transparency of the problem is preserved, as often the physics of the problem 
is lost among a great length of closed-form mathematical formulas; (2) future developments are 
simplified, as the basic computer program can be modified to describe different flow fields and/or 
different structures. 

A research program has therefore been undertaken to improve the analytical capability, and 
to develop practical computer programs for computing the displacement, velocity and stresses at 
different locations on realistic panel structures in different flow fields, including attached and 
separated flows at subsonic, transonic, and supersonic speeds. The first step in this process, 
discussed in this report, has been to consider the response of a clamped edge panel to the excitation 
of a subsonic attached turbulent boundary layer. The subsonic case was chosen because of the 
availability of corresponding excitation data (ref. 7) and response measurements (ref. 6). A similar 
analysis pertaining to attached and separated supersonic flow is in progress along with tests to 
obtain corresponding structural response data. For this reason, a computer program has been 
developed that is in modular form so that additions can be made to the program to accommodate 
different flow and structural conditions without affecting the rest of the program. The computer 
program pertaining to the analysis will be published separately. 


METHOD OF ANALYSIS 


Formulation of Displacement Power Spectral Density (DPSD) 

The displacement w(x,t) of a vibrating panel (fig. 1) is assumed to obey the classical thin 
plate equation: 

pw + cw + DV 4 W = p(x,t ) (1) 

In the present analysis a uniform plate is considered. Hence p and D are constants and independent 
of x. The above equation with c and p set equal to zero is the equation of free vibration of the 
plate, of which the solutions are called the natural normal modes of the plate. In the 

subsequent analysis, it will be assumed that the damping is so small that the natural normal modes 
and frequencies are not significantly changed. It is assumed further that xv(x,t ) can be expanded 
in terms of *p a (x) as follows: 
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( 2 ) 


w(x,t)=^q a (t)i/ a (x) 

Here \jj a is assumed to be properly normalized: 

J 'l' a (x)\l/p(x)dx = 5 a p 

A 

Furthermore, it will be assumed that \p a (x) can be obtained by separation of coordinates: 

^aC£) — n(x2 ) 

The term ip a (x) will be properly normalized if ) and ^ n {x 2 ) are normalized: 

I )d x i ~ &mn 

*2 

j i 1 )$n( x * )dx 2 — &mn 

J o 

It can be shown that, for a panel with simply supported edges 

FT mirxi 

*m( x i) = Jj: sm ~J~ i=1 > 2 

For a panel with clamped edges, \ l/ m takes on different forms according to whether m i 
or odd: 

(a) If m is odd, 


* m(Xi) = Tit if l cos ‘ ^ + km 008,1 Tm @ ' 3)1 

where y m are the roots of the equation 


. ym , . , ym „ 

tan — + tanh -y = 0 


and 


km ~ 




( 3 ) 
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(b) If m is even, 


$ nfai) 



+ k m sinh y m 



where y m are the roots of the equation 


(5) 


. 7 rn u ym 

tan — - tanh — 


= 0 


The orthogonality condition of the mode shape function can be used to find the normalization 
factors A m . 


A m ~ 


\^2y m 
Zm 


where 


Z m y?n ^ y^n k m 2 

~ y m ~~ s * n y iti + k m 2 
The numerical values of A m and y m are: 


(sinh y m ym) 

if m is odd 

(sinh ym — ym) 

if m is even 


m 

"4 m 

7m 

i 

0.7133 

4.730040 

2 

.7068 

7.853202 

3 

.7071 

10.995608 

4 

.7071 

14.137164 

5 

.7071 

17.278758 

6 

.7071 

20.420352 


Note that 


A m ~ 0.7071 


for m > 3 


and 


y m — y 6 + ( m - 6) 7 t for m > 1 

The generalized coordinate q a (t) satisfies the Lagrange equation 

m a q a {t) + C Q ^ a (0 + K a q a (t) = /<*(/) 
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where the generalized mass 


m a = n J \ l/ a 2 (x)dx = ju 


in view of the orthonormality condition of 4/ a (x). 
Similarly, 


C a = c 


Ka. — wiq [Uq ; 2 


and 


fa(0 = J P(x,t)\l/ a (x)dx (6) 

The above Lagrange equation is valid only for viscous damping. With structural damping the 
corresponding Lagrange equation (ref. 6) takes the form: 

m a q a (t) + K a (l + iv a ')q a (t ) = /«(/) 

where v a ' is the loss factor corresponding to structural damping. For harmonic solution the above 
equation can be cast into viscous damping form because q = ( l/icj)q . Therefore 

m a qa(0 + C'aeqait) + K a q a (t ) = f a (t) 

where C a e = {K a v a ’/u}) is the equivalent viscous damping coefficient. In practice, both viscous 
and structural damping are present, and the Lagrange equation takes the final form 

m a q a {t) + C ae q a (0 + K a q a (t) = f a (t ) (7) 


where 


Qtc — + C a 

is the equivalent damping coefficient that takes into account the effect of both viscous and 
structural damping. The viscous loss factor S a ' is usually defined by the following relation: 


Hence, 


_ Ka^oi j_ 

<-ae ~ T 


2K a S a ' 


co a 
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If the combined damping is small, so that the PSD curve is sharply peaked at the natural fre- 
quencies, only the frequencies around the neighborhood of the natural frequencies are of interest. 
Therefore 


CO ~ co a 


and 


Cae ~ 77 (25q/ 

CO 


W) 


CO 


where 


Voc ~ v a + 25a' 

is the combined loss factor for both structural and viscous damping. 1 It is this combined loss 
factor v a that is measured experimentally. The natural frequency co a will be computed by the 
method of Hearman (ref. 8). 

The Fourier transform of the Lagrange equation (7) takes the form: 

Q a (co) = /7 a (co)F a (co) (8) 


where 2 


H a (w) ma [( CJa 2 - cj 2 ) + iv a g ) a 2 ] 


(9) 


Equations (2) and (8) together give 

W(x,to) = HnicoWnicoH a(x) 

ft 

The displacement power spectral density is related to Wj{x,u) by 

Srfte.u) = ^lim j: Wt*(x,o))Wj(x,u) 


( 10 ) 


(ID 


'The loss factor v a takes into account the effect of hysteretic damping as well as viscous damping. The 
latter includes the interaction of panel motion on the flow field within and outside the boundary layer, and 
the effect of acoustical radiation into the interior of the structure (cavity), and into the exterior flow field. 
The present state of art requires that v a be measured experimentally. Development of analytical and experi- 
mental methods for the determination of loss factor are in progress. 

2 One can alternatively write Cae ^ where 25 q. = 26o( + r'a . This will give //(to) - 1 /wa [(wa - co ) 

+ 2i5aa>a<*>] instead of equation (9). The two expressions are equal at = co a . 
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where Wj is the truncated Fourier transform of w. Equations (7), (8), (10), and (11) together give 


Sjix.u) = E 5 )$ ol(x)4* 

a /S 

• eft’ dt" dx' dx" 


The excitation cross correlation function Rp(x',x",t ) and the cross spectral density Sp(x',x", o>) 
are defined by: 


Rp(x',x",T) 


1 , J 

lim ^ / p(x',t')p(x",t' + r)dt' 

7’->oo 


(12) 


S p (x’,x",a>)=j i - J°°R p (x',x",T)e- iuT dT 


(13) 


Equations (12) and (13) can be used to derive the following expression for Srf: 

S^(x,oj) - ASp(co) ^ ^a 2 (^)l^a( a? )l 2 /aa( w ) 
a 

+ 2AS p (oj ) X) + a (x)tQ(x)H a (u)H (j *(u)J a(} (u>) (14) 

a*j3 

a,/3 counted once 


where 


= ^ 5 ^" ' ) / / ta(xltp(x")Sp(x',x",u)dx' d£ 


A A 


(15) 


is the familiar structural acceptance first introduced by Powell (ref. 9). Note that it is dimensionless. 
An empirical formula for Sp(x',x", a>) is given in appendix A, and in appendix B the structural 
acceptance J a p is shown to be equal to the transition amplitude per unit area per unit excitation 
between the a and j3 mode of vibration. 


Separation Into Real and Imaginary Parts and Factorization of Acceptance 

Since Sp(x',x”,u ) is a complex function (see appendix A), the structural acceptance is 
generally complex. Therefore it can be written as 
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/a/3 i a (3 (3 


(16) 


Furthermore, with the assumption that both the mode shape function and the excitation can be 
factored into their respective coordinate components, the acceptances can be factored in the 
following forms: 


( q /3 inshnr 
kap = j' ns k m r 


for homogeneous turbulence 

1 f c > f Cl 


/ mr (co)=Y| f 'J 1 t m toc i , W r (x t ")\p p (li l ,0*>)\cosbi 1 dxi ' dxi 


i rS 2 

j' ns (u) =Y J J 'f / n(x2 , )'l's( x i , ')lf ) p( 0 ’Z2 dx i" 


-'o o 


1 r Cl 

Since the correlation coefficients Pp(£i ,0,u>) and pp(0,^ 2 ,w) do not depend on the signs of and 
£ 2 , that is, 

Pp(!.,0,«) = Pp(l*il.0,«) 

Pp(0,?2,w) = Pp(0,1^2 l.w) 

The above expressions for the acceptances can be simplified considerably if the region of integration 
is divided into two triangles bounded by the straight lines x” = 0, x' = £ and x = x ,x - 0, 
x" = £ and x" - x' and the symmetry property of the mode shape function is used: 

* m (*) = (-l) m+1 ^(£-*) 




The final results are: 



= 0 


if m + r even 


(17) 


if m + r odd 
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if n + s even 


£ x 

fnsM B Y 2 I 2 1 ‘ MOM*a")lV 0, * a,Cd)kfofj ' ' dx 2 " 
= 0 


if n + s odd 


( 18 ) 


kmr(u)=rj J * *m^i"^r(.Xi’ , )\Pp(iifi^)\sinbi l dx l , dx l n 

Inn 


0 0 


= 0 


if m + r odd \ 


if m+r even 


( 19 ) 


A physical argument leading to the second identity of the above equations will be given in 
appendix B. 


Two-Parameter Formulation of Acceptance 


As formulated above, the acceptances depend on the panel dimension C, , C 2 , the convection 
velocity U c , the frequency go, and the boundary-layer displacement thickness S*. Previous 
mathematical procedure requires recomputation of acceptances when any one of the above four 
parameters is changed. Since this recomputation is extremely time consuming, the acceptances 
have been reformulated as follows so that they do not depend on all of the above parameters. 
Define: 


Ft 


= Mi 

U c 


y t = 


Xj_ 

Ci 


71 i = y{ - yt 


/ / 


A 



where 


/= 1 , 2 


The acceptances then take the form: 
1 1 

i yi 


i jr l t _ p \ 

W (F *) = 2 J J T m (yi'Nr(yi")Pp 


0 0 
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cos nF 2 n 1 dyi dy x " ifw + reven^ 
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( 20 ) 


9 



dy 2 dy 2 ” 


if n + s even 




j' ns (F 2 ) = 2 J f ’)fs(y* ") P p {o,V 2 ?~Y ) 


o o 


= 0 


if n + s odd 


( 21 ) 


r r y ' - 

kmr(F \ ) ~ 2 J J i) m O'i ')>/vOT ") 
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sin 171 dy x ' dy\" ifm + /-odd 
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( 22 ) 


Here m {y) are the normalized mode shape functions (eqs. (4) and (5)) with £ replaced by unity; 
the pressure correlation coefficients \p p (vi, 0, itFJ 2)|, \p p (0, r? 2 , nF 2 /2)\ are as defined in 
appendix A. In equations (20) to (22), the acceptances now depend on only two parameters, 
Fi and A,-*. As long as A/* remains constant, the same acceptance curves can be used to compute 
the displacement power spectral density of different rectangular panels at different free-stream 
velocities. 


In terms of these new acceptances, the displacement power spectral density takes the form: 


S d (y,f) =S p (f)Y,* a 2 (y)\H a {2irf) P j mn {F l )f m (F 2 ) 


$y)his^2 )[£<* ) + hapkmr(Fi )] 

+ hap) (23) 

where 

g aP = -“> 2 )(ojp 2 -u 2 ) + v a i>poj a 2 iop 2 (24) 

= v a wo( 2 (w^ 2 - co 2 ) - v qu) p 2 (to Q 2 — to 2 ) (25) 

Acceptance charts based on equations (20) to (22) have been developed for the case of a 
rectangular panel under excitation of a turbulent boundary layer in a subsonic flow with a zero- 
longitudinal mean pressure gradient. The application of these charts in displacement PSD compu- 
tation is described in appendix C. 


+ 2 ^> 0 5 * 


RESULTS 


The first part of the results from the present analysis includes the response of a clamped edge 
panel to turbulent flow excitation applied on one side of the panel. The panel properties and flow 
parameters used for an illustrative computation are given in table 1. The results are shown in 
figures 2 to 9 along with some results of an identical panel with simply supported edges. Fig- 
ures 2 to 6 show the structural acceptances for a boundary layer with 5*/£ ( = 0.04468 and 
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= 0.06499, for which the corresponding response data are available from Wilby s experi 
ment (ref. 6). Figure 6 shows the variation of acceptance with boundary-layer displacement 
thickness. Figure 7 shows the calculated DPSD at the quarter- and midpoints of the panel, 
together with the corresponding experimental results. Figure 8 depicts the perspective view of 
the distribution of displacement PSD components at each mode on the panel. Figure 9(a) shows 
the variation of the displacement PSD with longitudinal distance x,/R, , evaluated at the lateral 
position x 2 /e 2 for which the PSD is a maximum for a given mode shown in figure 8. Similarly, 
figure 9(b) shows the variation of the displacement PSD with lateral distance x 2 /C 2 evaluated at 
the longitudinal position *,/£, for which the PSD is a maximum. PSD results are obtained by 
truncating the series in equation (23) after the seventh mode. It is found that decreasing the 
mode number to five, or increasing the mode number to nine, has a negligible effect on the results. 

The second part of the results (figs. 10 to 12) are charts of acceptances plotted against the 
frequency parameter Fj for various A /*. These charts, together with equation (23), enable the 
displacement PSD of the response of a clamped edge panel to be computed for attached turbulent 
boundary-layer excitations in subsonic flow. 


DISCUSSION 


Structural Acceptances 

Coincidence of structural and pressure waves - The response of a panel to a spatial-temporal 
correlated random pressure fluctuation caused by a turbulent boundary layer is characterized by 
wave-length matchings between the pressure wave and the flexural waves of the panel; a condition 
often referred to as “coincidence.” It is this peak coincidence that gives rise to large response if it 
occurs at one of the resonant frequencies of the structure. Coincidence frequencies are identified 
by peaks of the longitudinal joint acceptance curves (fig. 2). They occur at F (or 4/C/(/c) — 2m 
for both simply supported and clamped edge panels. 

Figure 2 shows that the degree. of wavelength matching is highest in the first mode, and 
decreases with increase of mode number for both boundary conditions. A simply supported panel 
exhibits a higher degree of matching than a clamped edge panel. Since the joint terms (the first 
sum on the right-hand side of eq. (23)) involving the joint acceptances account for almost all of 
the S^(x,co), it follows that a simply supported panel will undergo a larger mean square displace- 
ment than a clamped edge panel. 

In the case of a simply supported panel, the peaks occur at F < 2m because of the finiteness 
of the panel. This means that the flexural wavelength of the panel is shorter than the matching 
wavelength of the pressure wave. As the mode number is increased, the peaks approach F — 2m, 
the matching condition of an infinite panel, because the panel appears to be infinite when 
compared with the small wavelengths of the pressure waves at higher mode numbers. 

When the panel edges are clamped, the panel boundary has two effects, one due to the 
finiteness of the panel as stated above, and the other due to the rigidity of the clamped boundary, 
For a nonresponsive region of a clamped-edge panel the effect appears to be small when compared 
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with the long pressure waves at the lower modes, and large when compared with the short pres- 
sure waves at higher modes. Thus, the effect of the clamped edges is negligible at lower mode 
matchings. The acceptance curves (fig. 2) show that the peaks at lower mode numbers occur at 
F < 2m, because finiteness of the panel dominates. At higher mode numbers the opposite effect 
of the nonresponsive part of the panel edges begins to dominate and the peaks shift to F > 2m. 
It is significant that the peaks shift a half wavelength larger than 2m. Thus, the generally accepted 
assumption that wave matching always occurs at F = 2m is not true for a finite panel. 

Structural acceptances as transition amplitudes — Suppose an arbitrary external force is 
applied on a structure initially vibrating in the 0 normal mode. The structure will generally 
transit into another mode (the a mode) of vibration, which can be a normal mode or a linear 
superposition of normal modes of vibration. The probability of this transition is called the 
transition amplitude between the 0 normal mode, and the a mode when the structure is excited 
by the specific force in question. It is shown in appendix B that the structural acceptance J a p(co) 
is equal to the transition amplitude per unit area per unit excitation from the 0 to the a mode of 
vibration when the panel is excited by randomly fluctuating pressure forces. 

In the case of two-dimensional panel vibration the transition amplitude is regarded as a 
product of longitudinal and lateral acceptances. The longitudinal acceptance j mr - ik mr corresponds 
to the transition amplitude between the wth and the rth mode of vibration of a one— dimensional 
structure (a beam) under the excitation of the longitudinally correlated component of the pressure 
force or |p(£, , 0, co)|e 103 i u c i n equation (A 1). The in-phase part of this longitudinal component 
of the pressure force gives rise to the real part j mr of the longitudinal acceptance and the out-of- 
phase part to the imaginary part k mr . The lateral transition amplitude j' ns corresponds to the 
transition amplitude of another one-dimensional structure under the excitation of the laterally 
correlated component of the pressure force or lp(0,£ 2 ,w)|. 

Several interesting characteristics of the transition amplitudes are disclosed by the acceptances 
in figures 3, 4, and 5. Figures 3 and 4 show that longitudinal transition amplitudes /' „ and kvm, 
are typically maximum when m = n . In case m =£ n the magnitude decreases with increasing order 
of the sum of m and n. This implies that the probability of the initial mode remaining in the same 
mode is the highest and that there is less probability of transition to other modes. The figures also 
show that the acceptances have maxima at F— twice the lower mode number of m or n, indicating 
that the maximum transitions occur also at coincidence frequencies of the lower mode numbers 
of m or n. 

Figure 5 shows that the lateral acceptances f m are small compared with the longitudinal 
acceptances in figures 3 and 4, particularly at high mode numbers. Since the pressure field is 
less correlated in the lateral direction, it follows that the lateral acceptances should contribute 
less to the panel response than longitudinal acceptances (see eq. (23)). The appearance of wave 
matching in the lateral acceptances is not possible because of the lack of convection velocity in 
the lateral direction. 

Acceptance as a function of F and A*— Typical acceptances as shown in figure 6 vary signi- 
ficantly with A j * at low Strouhal number F, , and their dependence on A,* gradually disappears 
as F, increases. The reduction in dependency on A, * occurs because high values of F t correspond 
to high frequencies and short pressure wavelengths, and the boundary layer therefore appears to 
be infinitely thick in comparison with the short pressure waves. This factor was also taken into 
account when the empirical formulas for the correlation coefficients were derived as shown in 
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appendix A, where the correlation coefficients are independent of A,* at high frequencies. This 
same variation of wavelength with F, relative to A,* causes the acceptances for lower modes to 
vary more significantly with A , * than those involving higher modes. This variation can be seen 
by comparing typical acceptance j 2 2 with / 77 in figures 6(a) and 6(b), respectively. 


Structural Response 


Displacement PSD as function of frequency- Figures 7(a) and 7(b) show the results of the 
computed displacement power spectral density (DPSD) at the quarter point of the panel at 
M x = 0.3 and 0.5, and' their comparisons with Wilby’s measurements (ref. 6). Results are shown 
fo7 both clamped edge and simply supported panels. The dimensions and material properties of 
the panel are listed in table 1 . Figure 7(c) shows the results of the calculated panel response at the 
midpoint of the panel, where no experimental data are available for comparison. 

The study of figure 7 discloses that the peaks of the response curve always occur at, or very 
near, the natural frequencies. The envelope of the peaks in the spectra decays rapidly with fre- 
quency to justify the normal mode approach used in the present analysis. The response compu- 
tation based on the clamped edge boundary condition (figs. 7(a) and 7(b)) agrees better with 
Wilby’s measurements (ref. 6) than with the computations based on the simply .supported boundary 
condition. The results also show that the joint terms (the first sum of eq. (23)) account for almost 
aU of the contribution to DPSD at the peaks. The contribution of the cross terms (the second 
sum of eq. (23)) to DPSD is completely negligible except at some valleys of the DPSD curve. 
Powell (ref. 9) showed that because of the orthogonality condition of the mode shape functions, 
the cross-term contribution to a uniformly loaded surface will be zero when averaged over the 
entire surface. However, he also pointed out that the cross terms cannot be neglected if damping 
is appreciable or if the force is applied to a localized area. 

When the Mach number is changed from 0.3 to 0.5 (figs. 7(a) and 7(b), the overall shape of 
the DPSD curve does not change, but the displacement per unit excitation increases at higher 
frequencies. The corresponding rms displacement per unit pressure input increases by about 
11 percent at the quarter point and 17 percent at the midpoint of the panel. This increase in 
response is attributed to the increase in the spatial correlation of the pressure field at high fre- 
quencies as the Mach number is increased. Since pressure excitation also increased with velocity 
( N /p - =« 0.006 q^,), the rms displacement of the panel will increase, as a result of both the increase 
in pressure excitation Sp(f) and the increase in DPSD per unit pressure excitation, S^(x,f)/Sp(f). 

The computed DPSD at the midpoint of the panel (fig. 7(c)) does not have peaks at natural 
frequencies, which correspond to antisymmetric modes, because the panel midpoint is a node 
point for the antisymmetric modes. 

Displacement PSD distribution on a panel- Computer programs capable of displaying the 
perspective view of displacement PSD on a panel were also developed. The displacement PSD on 
a panel is necessary information for the assessment of the distribution of stress. Typical results 
given in figure 8 show the distribution of PSD components on a clamped edge panel at natural 
frequencies corresponding to modes (1-1), (2-1), (1-2), (3-1), and (2-2). Figure 9 shows the 
variation of maximum displacement PSD with longitudinal and lateral coordinates. Note that 
when mode m or n > 1 , the peaks of the PSD next to the panel edge do not occur at exactly 
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i/2m from the edge, as they would in the case of a simply supported panel. The fact that the 
peaks of the PSD occur at distances greater than 8/2 m from the edge is due to the rigidity of the 
panel at the clamped edge. Note also that when mode m or n > I and the mode number is odd 
(for instance m = 3 in fig. 9(a)), the PSD in the middle region of the panel is symmetric, but the 
PSD next to the panel edge is asymmetric about the respective peaks. 


CONCLUSIONS 


The study on the response of a rectangular panel under the excitation of a turbulent boundary 
layer discloses the following: 

1. The matching between the pressure wave and the flexural wave of a panel does not 
occur exactly at even integer values of the reduced frequency F = 4/8/(/ c , as it would when the 
panel is infinite for two reasons; (a) the finiteness of the panel significantly affects the low fre- 
quency matchings, and (b) the rigidity of the clamped boundary significantly affects the high 
frequency matchings. 

2. The structural acceptance J a p is equal to the probability of transition, or transition 
amplitude, from the 0 mode to a mode of vibration, per unit excitation power spectral density 
per unit area of the panel, when the panel is excited by the random pressure. The longitudinal 
transition amplitudes (j mn and k mn ) exhibit a general behavior, being maximal when m = n. 
When m =£ n, the magnitude decreases with increasing order of the sum of m and n. This decrease 
implies that the probability of preserving the same mode is higher than the probability of transiting 
to different modes. 

3. An expression that depends only on the dimensionless frequency F/ and boundary 
layer thickness A/* can be derived for the structural acceptance. Corresponding acceptance 
charts were therefore prepared which make possible a quick estimation of the response of any 
rectangular panel to excitation of a turbulent boundary layer in subsonic flow with a zero 
longitudinal mean gradient. 

4. The lateral acceptances are small compared with the longitudinal acceptances. Since 
pressure field is less correlated in the lateral direction, the lateral acceptances contribute less to 
the panel response than longitudinal acceptances. 

5. The response computation based on the clamped edge boundary condition of a panel 
is in better agreement with the measurements than previous results obtained by investigators who 
used a simply supported boundary condition in the acceptance calculation. 

6. As the Mach number is increased within the subsonic range, the response of a panel is 
also increased in its high frequency components of vibration. The increase is attributed to two 
factors: the increase in the spatial correlation of the pressure field at high frequencies; and the 
increase of the magnitude of pressure excitations. 
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7. The variation of boundary-layer thickness affects only the low frequency component 
of the response through the structural acceptances, which vary significantly near the low fre- 
quency, particularly where the wave matching occurs. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, July 7, 1972 
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APPENDIX A 


RANDOM EXCITATION BY SURFACE PRESSURE FLUCTUATIONS 


The pressure fluctuations that occur on a surface adjacent to a turbulent boundary layer 
have been shown by many investigators (refs. 7 and 10 to 16) to be spatially and temporarily 
correlated. For homogeneous turbulence, as attached flow, the cross spectral density of the 
fluctuating pressure depends only on the separation distance and can be expressed in the following 
form: 




Sp(& i &2 A>) 
Sp(oj) 

= 1 (Al) 

where Ip^d,, £ 2( cj)I is the cross correlation coefficient between two points separated by distances 
of and £ 2 in the .*! and x 2 directions. Experimental data also show that the correlation 

coefficient can be separated in coordinates as follows: 



IPp(tl 2 A>)| = 

|p p a J ,0,co)| lpp(0,| 2 , CO )| (A2) 

In the case of a turbulent boundary layer at a zero longitudinal pressure gradient in a subsonic 
flow, empirical formulas have been constructed for |pp(£, ,0,o;)| and \pp(0£ 2 ,ta)| based on Bull’s 
experimental data (ref. 7): 


Pp(vi Ay 1 ) 

7 j p . 

s e -a t y 1 It? i 1 

if yp > ki 



= e -0Ci \V i I/a I * 

ify 1 A,* < A;, 


Pp(o,r?2 .y 2 ) 

j tF 2 . , 

= e ~^3~ |T?2 1 

if |t? 2 I > - [9.1 log (^y A 2 *) + 5.45] A 2 * 



= c + de~°* 1/A 2 * 

if It ? 2 1 < - [9.1 log (y 1 A 2 *) + 5.45] A 2 * 

where 





a, = 0.1 

c = 0.28 



a 2 = 0.037 

d = 0.72 



a 3 = 0.715 

k t = 0.37 



a 4 = 0.547 
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The convection velocity, U c , in this mathematical model is given by, 

U c = (k , + k 2 <T a swS*/f/oo )f/oo ( A3 ) 

where 

k i = 0.59 
k 2 = 0.30 
a 5 = 0.89 
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APPENDIX B 


ACCEPTANCE AS TRANSITION AMPLITUDES BETWEEN TWO NORMAL MODES 
OF VIBRATION AND ITS SIGNIFICANCE IN PANEL VIBRATION 


In view of the importance of acceptances in predicting structural response to random forces, 
a physical meaning of acceptances is derived below. Suppose a structure is vibrating initially in 
the 0 normal mode. If an arbitrary external force is applied, the structure will in general transit 
into another mode of vibration which can be a linear superposition of normal modes or a pure 
normal mode. The probability of this transition is called the transition amplitude between the 
(3 normal mode and the final mode when the structure is excited by the specific force in question. 
Since the final mode is, in general, a linear superposition of the normal modes, it can be expressed 
as S a aP'l / oi> where a a $ is the transition amplitude between the a and )3 normal mode. 

To show the relationship between the transition amplitude and the acceptance, a mathematical 
formulation of the above discussion is necessary. In the language of linear algebra, the method of 
normal mode analysis assumes that the normalized mode shape functions \p a form an orthonormal 
basis in a linear space, S, which is the space of all possible vibration modes of the structure. Any 
mode of vibration is a vector in S and therefore can be expressed as a linear sum of i p a . The 
orthonormal condition of \ p a can be expressed mathematically as: 

J A 'P aQD'P p(x)dx = 8 a p (Bl) 


where the integration is over the entire structure, or in the present case, the area of the rectangular 
panel. This suggests that the scalar product between two vectors A, B in S can be defined as: 


(A\B) = 


l 


A {x)B(x)dx 


(B2) 


In this linear space, a disturbance such as Sp(x',x",w )is described mathematically by a linear 
operator 0£^(x',x'', co) in S with the property tnat if A is a vector in S, 


O s (x',x",io)A = f dx" S p (x',x",cj)A(x") 


(B3) 


in accordance with the definition of scalar product in equation (B2). By definition of O $ , O $ A 

is also in S. In particular, when A = O s \p p is also in S. Therefore it can be expanded^as 

follows: P 


In the case of a flat panel, each of a,(3,y contains two indices such as (1,1), (1,2) . . . ( oo , o °).^The 
coefficients a a(3 can be obtained by forming the scalar product between \p a and 0 Sp (x',x",u ) 
as follows, 


W a(x')\Os p (x\x'' p(x'')) = <'Pa(x')\^,a y ^\p J (x^)) 

= E a 70 <\M*')I^ 7 (*')> 




ay 


a a($ 


By definition of equations (B2) and (B3), 

a(x r )\Og n (x\x'\G})\p p(£ r )) = f [ 4* a(x')Sp(x\x’ T )dx dx_ 

A A 


(B4) 


(B5) 


Comparing equations (B4) and (B5) with the definition of the acceptance 


j a p( co) = s (lyA jf^cc(xlS p (x\x'\u)4 /e (x'')dx' dx" 


the following identification can be established: 


/ a 0(w) Sp ^ A a ap 


Thus the acceptance J a p is just the transition amplitude per unit PSD of excitation per unit area 
from the j3 to the a normal mode of vibration when the structure is excited by the disturbance 
described by Sp(x',x",o>). 

In the case of two-dimensional panel vibration the transition amplitude is a product of 
longitudinal and lateral acceptances (see equations (14) and (23)). The longitudinal acceptance 
j mr -ik mr corresponds to the transition amplitude between the m and the r mode of vibration 
of a one-dimensional structure under the excitation of the longitudinal correlated components 
of the pressure force. The amplitude of a final mode r formed from initial modes 1, 2, ... n is 
proportional to the following sum: 


n 


E 

m = l 


(j mr + kmr) 


where the first sum represents the transitions caused by the in-phase pressure fluctuations, and 
the second sum by the out-of-phase pressure fluctuations. When any two of the above acceptances 
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have the same sign at a given frequency, the transition is such that they reinforce each other, 
giving rise to constructive interference. If they have opposite signs, they cancel each other, giving 
rise to destructive interference. 

During a transition, the parity of the mode must be conserved if the excitation has even 
parity or symmetry in space, and must change if the excitation has odd parity or antisymmetry 
in space. Since j mr is proportional to transition amplitude caused by the in-phase pressure force, 
which is symmetric in space (it contains a cosine term in eq. (20)), j mr must be zero if m + r is 
odd. On the other hand, since k mr is proportional to the transition amplitude caused by the 
out-of-phase pressure force, which is antisymmetric in space (it contains a sine term in eq. (22)), 
k mr must be zero if m + r is even. For the same reason f mr must be zero if m + r is odd. The 
physical significance of equations (20) to (22) is thus shown in terms of the parity of modal 
transition. 
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APPENDIX C 


APPLICATION OF ACCEPTANCE CHARTS IN 
DISPLACEMENT PSD COMPUTATION 


Values of j mr {F x ), k mr (F i ), and f ns (F 2 ) up to the seventh mode have been computed and 
plotted in figures 10 through 12 for F up to 15.0, and A/* ranging from 0.015 to 0.08 in incre- 
ments of 0.01, for a clamped edge rectangular panel under boundary layer excitations. These 
charts, together with equation (23), enable the displacement PSD be computed for a clamped 
edge rectangular panel exposed to a subsonic freestream velocity. The procedure is to compute 
the U c first, using equation (A3). Then, knowing the dimensions of the panel and the boundary 
layer displacement thickness, the two parameters F[ and A/* can be computed. The corresponding 
acceptances can be read from the charts and substituted into equation (23) with natural frequencies 
and loss factors which can be determined experimentally or analytically. The method of linear 
interpolation can be used if A/* does not coincide exactly with one of the values at which the 
acceptances are evaluated. 
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TABLE 1 PANEL PROPERTIES AND FLOW PARAMETERS 

Ref. 6 

Mild steel panel 

Cj 0.1016 m (4 in.) 

£ 2 0.06985 m (2.75 in.) 

h 0.000381 m (0.015 in.) 

E 0. 2323X10* 2 N/m 2 (33.7X10 6 lb/in. 2 ) 
o 0.3 
v 0.009 

Density 

7473 kg/m 3 (0.27 lb/in. 3 ) 

Flow parameters 

Moo 0.3 and 0.5 
5* 0.00454 m (0.179 in.) 




Figure 2.— Longitudinal joint acceptances of clamped and simply supported panel for A* = 0.04468 
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Figure 3. — Continued. 
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Figure 7. — Displacement PSD of the panel response to unit excitation PSD. 
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Figure 8.- Perspective view of displacement PSD of the response of clamped edge panel to unit excitation at M = 0.3. 
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Figure 8. - Continued. 






Natural frequency 
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Figure 8. — Continued. 
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Figure 8. - Concluded. 
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(b) Variation with lateral coordinates. 
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Figure 10. — Real part of longitudinal acceptance of clamped edge panels. 
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Figure 10. - Continued 
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Figure 10. — Continued. 
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Figure 10. - Continued. 
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Figure 10. — Continued 









Figure 10. — Continued 
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Figure 10. — Concluded. 
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Figure 11. — Continued. 
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Figure 1 1. - Continued. 
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Figure 11. — Continued. 
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Figure 11.— Continued. 
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Figure 11. — Concluded. 
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Figure 12. — Continued. 
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Figure 12. - Continued. 




Figure 12. — Continued. 
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Figure 12. - Continued. 
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Figure 12. - Continued. 
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Figure 12. - Continued. 
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